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Abstract. From simulations of transit observations, it is found that the 
detectabihty of extrasolar planets depends only on two parameters: The 
signal-to- noise ratio during a transit, and the number of data points ob- 
served during transits. All other physical parameters describing transit 
configurations (planet and star size, orbital period, orbital half axis, lat- 
itude of the transit across star) can be reduced to these two parameters. 
In turn, once the dependency between transit detectabihty and these two 
parameters has been determined, predictions for an instrument's ability 
to detect transits by any combination of physical parameters can be de- 
rived with ease. These predictions are applied to the Eddington proposal 
of a combined Astroseismology/Transit-detection space mission currently 
under study by ESA, which is described briefly. 



Introduction 



To evaluate the performance for the detection of extrasolar planets of observa- 
tions with the transit method, Jenkins, Doyle & Culler (1996) introduced the 
method of the calculation of two distributions, one of the detection statistics for 
transit-less noise-only data (which may be observed or simulated), and one de- 
tection statistics for data with known transits in them. This method allows the 
determination of detection probabilities while at the same time setting limits for 
the false alarm rate - important in the case of transit detections, as one neither 
knows if a transit signal is in the data, nor the signal's duration, amplitude or 
period. Doyle et al. (2000) used this method to set detection limits in ground 
based observational data, taken with the goal to detect transits across the eclips- 
ing binary CM Draconis. The evaluations shown here were performed during 
the course of specifying the Eddington proposal currently pending at ESA for 
a combined astroseismology /planet-finding mission with a Im-class space tele- 
scope, which is described at the end. In the following, the dependence of the 
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parameters describing a transit configuration (Planet's size, star's size, planet's 
period for example) on detectability is shown. 

2. Extrasolar planet detectability by transits 

Planet detection performance may be expressed in terms of detectability or the 
detection probability pd, which is the probability that a physically present tran- 
siting planet is going to be detected by the experiment. The detectability of 
a planet in an observed transit sequence depends on a number of physical pa- 
rameters of the planet-star system and of the observing instrument. The most 
obvious one is the size of the planet, with the amount of brightness-loss during 
a transit being proportional to the square of the planet radius, Rpi. Similar, 
detectability will remain roughly constant, if the ratio Rpi/R^, is kept constant. 

To asses the performance in the detection of planetary transits, simulated 
lightcurves with the instrument's noise characteristics were submitted to a de- 
tection test. This detection test gives some 'detection value', which may be 
calculated by any kind of procedure that performs some evaluation between a 
given model lightcurve and the data. The most suitable one for data which are 
dominated by 'white' noises and without (or with correctable) zero-point drifts 
is a scalar multiplication between the data vector and the model vector: 



where C is the detection value, and di, rrii are data and model values. 
The detection test was performed in large numbers on lightcurves with - and 
without simulated transits added to the noise The highest 'detection values' 
from a large number (between 10^ — 10^) of lightcurves without transits (that 
is, with noise only) give a threshold for detection values that constitute false 
alarms. The fraction of lightcurves with transits (again using a large number 
of randomly generated curves, but with known transit parameters) that results 
in detection values higher then the false alarm threshold gives the detection 
probability. This procedure is similar to one used for an assessment of detection 
limits in the ground based transit detection project 'TEP' (Doyle, et al. 2000). 

During the course of evaluating a variety of transit configurations, it be- 
came obvious that the physical parameters governing transit lightcurves may 
be reduced to two parameters which solely govern the detectability in any ob- 
served sequence of transits around single stars (the case for binaries would be 
more complicated): the data's S/N during a transit, and the total number of 
data points within transits, Nfp. Thus, if the function p^^S/N ,Ntp) can be es- 
tablished, the detectability can be calculated for any combination of physical 
parameters. The physical parameters influence S/N and Ntp as follows: 

- The Signal S depends on planet's and star's sizes, the latitude of transit, 
and on the star's limb-darkening (mostly for high-latitude transits). For 
a precise calculation of model transit lightcurves (and hence S), a code 
dealing correctly with the planet's ingress and egress and with differential 
limb-darkening under the planet's disk should be used, such as 'utm' (Deeg, 




1999). 



Detection limits for transit observations 



3 



- The noise N depends on the noises affecting the observing instrument's 
performance (which has always a dependency on a star's apparent bright- 
ness) and on intrinsic brightness fluctuations (stellar activity, oscillations, 
variability) 

- The number of on-transit data points, Ntp, depends on exposure-time, the 
number of transits observed (=observational coverage) and on the duration 
of a transit. The duration in turn is governed by the planet's orbital period, 
by the star's mass and size, by the latitude of the transit, and to a small 
amount by the size of the planet. It can be calculated by straightforward 
application of basic geometry and Kepler's third law. 

For detectability evaluations, a 'standard case' was chosen, which has a 
detection probability {pd) of 0.99. In Fig. 1, the sequence of 6 transits of the 
standard case is shown. That case comprised the following setup: 

- Physical parameters: KO star (0.85 Rsoi) 0.78 Mg^i), orbited by a planet 
with 1.4 REaxth) 200d orbital period, transiting at a latitude of 45°across 
the star. 

- Observational parameters: 3 years (1095 days) observing duration; lightcurve 
with a noise of 450 ppm, or approximately O.Smmag from 600 seconds ex- 
posure time. In the case of the Eddington proposal, this is expected for a 
star with 16 mag. 

As an example of parameter variations around the standard case, the dependency 
of Pd on planet-size is shown: 



Rpl /REarth 


S/N 


Pd 


1.0 


0.25 


0.16 


1.3 


0.42 
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1.4 


0.48 


0.99 


;^i-5 


^0.55 


1.00 



More useful then detection probability is however the determination of the 
minimum planet size that can be detected reliably. A 'reliable detection' is con- 
sidered to have pd > 0.99. From the table above, the detection limit is therefore 
a planet of 1.4 REarth- Similar tables where also generated for the 'standard 
case' but with orbital periods of 50, 100 and 365 days, thus performing an em- 
piric determination of the planet-size -Rpi,min where pd ~ 0.99. These results 
are shown in the next table, where besides -Rpi,min the following quantities are 
given: utr, the number of transits in the 3-year long lightcurve; tfr, the duration 
of one transit; Ntp, the total number of on-transit data points at 600 seconds 
exposure time; and S/N, the signal to noise ratio, where the signal is taken as 
the maximum relative brightness variation during a transit. 
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Figure 1. Left: Sequence of 6 transits corresponding to observations 
of the 'standard case'. The upper panel gives noiseless model transits, 
the lower one is the simulated noisy lightcurve with the model sub- 
tracted. Of the whole lightcurve, only segments near the 6 transits are 
shown. As the detection limit of pj=0.99 corresponds here to a S/N of 
0.5, individual transits arc not detectable - neither to the eye nor to any 
algorithm. Right: the same transit sequence, but data are rebinned by 
2 and all 6 transits are coadded (that means: the model-transit is 6 
times deeper, and the noise is ^6/2 times larger). Center panel is pure 
noise, lower panel is noise+transit. Planet detectability is practically 
independent of exposure time or binning of the data. This holds, as 
long as noises dominate that are proportional to (exp.time)~^/^, such 
as photon and background noise. 



From above table, a linear relationship between logs'/A^ and log Ntp was 
established, which can be expressed as: 

log S/N = 0.90 - 0.50 log Ntp, 

or more generally: 

S/N * ^J~{Ntp) = const for p^ = const. 

It should be noted, that this expression is valid only when using C = J2di * mi 
for the determination of the 'detection value'. Other possibilities to derive de- 
tection vahics, such as the calculation of values between model and data, are 
expected to lead to other dependencies, and may not allow the reduction of the 
physical transit parameters to the two parameters S/N and Nfp. Above expres- 
sion between the fundamental parameters governing transit detectability (using 
the level pd = 0.99) allows the determination of i?pi,min for any combination of 
physical parameters constituting a transit configuration. 

In Fig. 2, i?pi,min is given for a variety of stellar types and orbital periods 
and apparent magnitudes as expected by Eddington. As the calculations were 
performed for transits at latitudes of 45°, the given planet sizes are valid for 
70% of theoretically possible transit configurations. For the remaining 30% 
of polar transits, detection probabilities may be significantly lower, and planets 
with -Rpi^min cannot be detected reliably. Indicated as hatched regions in Fig. 2 
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Figure 2. Minimum detectable planet size (in -Rsarth) for a range of 
stars in dependence of orbital period (in days) and apparent magnitude 
for the Eddington proposal. A transit latitude of 45°was used. Hatched 
regions indicate the circumstellar habitable zone. 
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is the extent of the circumstehar habitable zone based on values from Whitmire 
&: Reynolds (1996). It can be seen, that Earth-sized planets around a solar- type 
star of rriy = 14 may be detected, and even stars with niy = 18 may still allow 
for the detection of 'large terrestrial' planets. 

3. The Eddington proposal 

This work was performed in the context of the Eddington study. Eddington is a a 
combined asteroseismology and planet-detection mission, based on a a 1 m class 
space telescope, initiated by a European consortium of scientists (Favata, Rox- 
burgh, & Christensen-Dalsgaard 2000) in response to ESA's F2/F3 mission call. 
Eddington was successfully selected for the initial currently ongoing study phase, 
and will undergo a final selection in September 2000. With foreseen photometric 
precision of ~ 1 part in 100 000, the major science goal for the planet-detection 
part will be the search for terrestrial size planets in the habitable zone. With its 
large mirror size (1.2 m diameter) and field of view (3° diameter), it will monitor 
several hundred thousand stars during the 3 years of its mission dedicated to 
the transit search. To enable long, uninterrupted observations within a stable 
thermal environment, an L2 halo orbit is being baselined. These parameters will 
allow Eddington to achieve qualitatively different science goals than the small 
missions currently in preparation (COROT, MONS, MOST). The discovery of 
extrasolar terrestrial planets, in addition to being a major achievement on its 
own, will serve as a fundamental stepping stone for future missions specifically 
designed for the detection of extraterrestrial habitats, such as ESA and NASA 
missions IRSI-Darwin and Terrestrial Planet Finder (TPF). Eddington will also 
contribute in a major way to the field of astroseismology, enabling an accu- 
rate characterization of the stellar structure for stars spanning a wide range of 
masses, ages and chemical compositions. 
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